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>. We perform a comparative study of the neutralino dark matter scattering on nucleon in three 


popular supersymmetric models: the minimal (MSSM), the next-to-minimal (NMSSM) and the 
nearly minimal (nMSSM). First, we give the predictions of the elastic cross section by scanning 
over the parameter space allowed by various direct and indirect constraints, which are from the 
measurement of the cosmic dark matter relic density, the collider search for Higgs boson and sparti- 
cles, the precision electroweak measurements and the muon anomalous magnetic moment. Then we 
demonstrate the property of the allowed parameter space with/without the new limits from CDMS 
II. We obtain the following observations: (i) For each model the new CDMS limits can exclude a 
large part of the parameter space allowed by current collider constraints; (ii) The property of the 
allowed parameter space is similar for MSSM and NMSSM, but quite different for nMSSM; (iii) The 
future SuperCDMS can cover most part of the allowed parameter space for each model. 
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I. INTRODUCTION 


Although there are many theoretical or aesthetical arguments for the necessity of TeV-scale new physics, the most 
convincing evidence is from the WMAP (Wilkinson Microwave Anisotropy Probe) observation of the cosmic cold 
dark matter, which naturally indicate the existence of WIMPs (Weakly Interacting Massive Particle) beyond the 
prediction of the Standard Model (SM). By contrast, the neutrino oscillation may rather imply trivial new physics 
(plainly adding right-handed neutrinos to the SM) or new physics at some very high see-saw scale unaccessible to any 
foreseeable colliders. Therefore, the TeV-scale new physics to be unravelled at the Large Hadron Collider (LHC) is 
most likely related to the WIMP dark matter. 

If WIMP dark matter is chosen by nature, then it will naturally direct to low-energy supersymmetry (SUSY) 
with R-parity although other miscellaneous speculations are also possible. In addition to the perfect explanation 
of cosmic dark matter, to make perfection still more perfect, SUSY can also solve another plausible puzzle, namely 
the 30 deviation of the muon anomalous magnetic moment from the SM prediction. In the framework of SUSY, 
the most intensively studied model is the minimal supersymmetric standard model (MSSM) {a}, which is the most 
economical realization of SUSY. Since this model suffers from the -problem and the little hierarchy problem, other 
supersymmetric models have recently attracted much attention, among which is the extension by introducing a 
gauge singlet superfield S, such as the next-to-minimal supersymmetric model (NMSSM) (a and the nearly minimal 
supersymmetric model (nMSSM) [3h W. In addition to the attractive phenomenological virtues like the alleviation 
of the little hierarchy problem and the possible explanation [s] of PAMELA positron excess (albeit subject to large 
uncertainty and could be explained astrophysically by pulsars) id, such singlet extensions are arguably motivated by 
some fancy string theory, e.g., the NMSSM can be constructed from a heterotic string [7]. 

In this work, motivated by the CDMS II new results E! id, we examine the SUSY dark matter scattering on the 
nucleon (x-nucleon scattering). In the literature such a topic has been studied mainly in the constrained MSSM ig. 
Our work is projected to have the following features: 


(i) We perform a comparative study for three popular SUSY models: the MSSM, the NMSSM and the nMSSM. 


(ii) We consider the constraints from the cosmic dark matter relic density and current collider experiments, such 
as the collider search for Higgs boson and sparticles, the precision electroweak measurements and the muon 
anomalous magnetic moment. By scanning over the parameter space subject to these constraints, for each 
model we find out the allowed parameter space and give the predictions of the cross section for y-nucleon 
scattering with comparison to the CDMS II results. 


(iii) We demonstrate the properties of the allowed parameter space (such as the components of the neutralino dark 
matter and the invisible Higgs boson decay into a pair of dark matter particles) by comparing the three models. 


(iv) We show the capability of the SuperCDMS in probing the currently allowed parameter space for each model. 


This paper is organized as follows. In Sec.II we briefly describe the three models: the MSSM, the NMSSM and the 
nMSSM, focusing on the Higgs sector and the neutralino/chargino sector since they are directly relevant to the dark 
matter scattering. In Sec.III we scan over the parameter space under current constraints, and give the predictions 
of the cross section for y-nucleon scattering with comparison to the CDMS II results. Also we will demonstrate the 
properties of the allowed parameter space with/without considering the CDMS new limits. In Sec. IV we give our 
conclusions. 


II. SUPERSYMMETRIC MODELS 


As the economical realizations of supersymmetry, the MSSM has the minimal content of particles, while the NMSSM 
and nMSSM extend the MSSM by only adding one singlet Higgs superfield S. The difference between these models 
is reflected in their superpotential: 


Wussm = Wr + pfu- Ha, (1) 
a e Ea 

Wnmssm = Wr + àHu : H4S + 35", (2) 

Wamssm = Wr +AHy- HaS + Er M29, (3) 


where Wp = YQ . HU — YaQ . HyD — Y.-L . Hab with Q, U and D being the squark superfields, and Ê and Ê being 
the slepton superfields, H,, and Hg are the Higgs doublet superfields, À, « and £p are dimensionless coefficients, and 
u and Mn are parameters with mass dimension. Note that there is no explicit -term in the NMSSM or nMSSM, and 
an effective u-parameter can be generated when the scalar component (S) of S develops a vev (vacuum expectation 
value). Also note that the nMSSM differs from the NMSSM in the last term with the trilinear singlet term KS? of 
the NMSSM replaced by the tadpole term € pM2S . As pointed out in Ep such a tadpole term can be generated at a 
high loop level and naturally be of the SUSY breaking scale. The advantage of such replacement is the nMSSM has 
no discrete symmetry and thus free of the domain wall problem which the NMSSM suffers from. 

Corresponding to the superpotential, the Higgs soft terms in the scalar potentials are also different for three models 
(the soft terms for gauginos and sfermions are the same and not listed here) 


Vote = ma|Hal? + m2|Hu|? + (Buby -Ha + hic.) (4) 
VNMSSM 2) Hal? + m2 | H,,|? + MISI + (ASH Hy + AnS? + he.) (5) 
VAROM = MAJHA + m2 ||? + m2|S)? + (AAS Ha: Hy + Es MBS + hic.) . (6) 


After the scalar fields H,,,Hqg and S develop their vevs v,,, vg and s respectively, they can be expanded as 


Z (va + ba + ipa) Hy 1 , 
m=(7 n d Eae so a ek) (7) 


The mass eigenstates can be obtained by unitary rotations 


Hı Qa Ay Pa 4 i 
Hə = Qu ; Ag = U^ Pu , ( ) = ge ( nd ) , (8) 
Az oO Go E M 


where Hi 2,3 and Aj. are respectively the CP-even and CP-odd neutral Higgs bosons, G? and Gt are Goldstone 
bosons, and H* is the charged Higgs boson. So in the NMSSM and nMSSM, there exist a pair of charged Higgs 
bosons, three CP-even and two CP-odd neutral Higgs bosons. In the MSSM, due to the absence of S, we only have 
two CP-even and one CP-odd neutral Higgs bosons in addition to a pair of charged Higgs bosons. 

The MSSM predict four neutralinos x? (i = 1,2,3,4), i.e. the mixture of neutral gauginos (bino à’ and neutral 
wino À?) and neutral Higgsinos (Y, , Yi, ), while the NMSSM and nMSSM predict one more neutralino because the 


singlino 7g comes into the mixing. In the basis (—id’,—iA°, Y% , WY, Ws) (for MSSM ys is absent) the neutralino 
mass matrix is given by 


Mı 0 mzswśsb >= mMzSw Cb 

i Ma Mamay. Mare for MSSM (9) 
MZSWSbh —MNZSWSb 0 — 
—MZSWwWCh —™MNZCWCH =u 0 

Mı 0 mzswśsb —MZSWwCh 0 

0 Mo —mzCcwsSb MzcwCp 0 
MzSwsSh —MzSw Sd 0 —u — ÀAUCp for NMSSM (10) 
—Mzswth —Mzcw Cp —u 0 —vSsp 

0 0 —vcp —vsp 254 

Mı 0 MZSWSh —MNZSWwCp 0 

0 Mo —Mzcwsp MzcwCp 0 
MzswsSh —MzswSp 0 —u —Avcy for nMSSM (11) 
—MzswCch —Mzcw Cp —p 0 —vsp 

0 0 —Xvcp —vsp 0 


where Mı and Mo are respectively U(1) and SU(2) gaugino masses, sw = sinOw, cw = cosOw, sẹ = sin’ and 
c = cos 6 with tan 8 = vu/va. In our study the lightest neutralino y? is assumed to be the lightest supersymmetric 
particle (LSP), serving as the SUSY dark matter particle. It is composed by 


x? = Nir(—tA1) + Nio(—iA2) + Miah, + Nady, + Nisvs, (12) 


where N is the unitary matrix (Nj5 is zero for the MSSM) to diagonalize the mass matrix in Eqs. (911). 
The chargino sector of these three models is the same except that for the NMSSM/nMSSM the parameter pi is 


replaced by Hef. The charginos x75 (m+ < myx) are the mixture of charged Higgsinos YH, P and winos AF = 


A! + à2)/ V2, whose mass matrix in the basis of (—iA+, y is given by 
Hua 


Mo V2mw sin 8 (13) 
V2mw cos 8 [left 


So the chargino x7 can be wino-dominant (when Mə is much smaller than u) or Higgsino-dominant (when u is much 
smaller than M2). Since the composing property (wino-like, bino-like, Higgsino-like or singlino-like) of the LSP and 


the chargino yj is very important for SUSY phenomenology, we will show such a property in our following study. 


Ill. NUMERICAL RESULTS AND DISCUSSIONS 


So far there are various constraints from both collider and dark matter experiments. In our study we consider the 
following constraints: 


(1) Direct bounds on sparticle and Higgs masses from LEP and Tevatron experiments ud, egm + > 103.5 GeV, 
me > 73 GeV, mg > 94 GeV, mz > 81.9 GeV and my+ > 78.6 GeV. 


(2) LEP II search for Higgs boson a3), which include various channels of Higgs boson productions {14}. 


(3) LEP I and LEP II constraints on the productions of neutralinos and charginos, including the LEP I invisible 
Z- mel T(Z > x?x?) < 1.76 MeV, the LEP II neutralino production o(ete~ + y$x?) < 107? pb (i > 1) and 
a(ete™ — x?x}) < 107} pb. 


(4) Indirect constraints from precision electroweak observables such as pe, sin” of rp and My, or their combinations 
cili = 1,2,3) [15]. We require e; to be compatible with the LEP/SLD data at 95% confidence level. Also, for 
Ry =T(Z > bb)/T(Z — hadrons) whose measured value is Ry? = 0.21629 + 0.00066 and the SM predicion 
is RIM = = 0.21578 for m = 173 GeV E we require RUSY is within the 20 range of its experimental value. 
Ver i B-physics constraints are also included 


(5) Indirect constraint from the muon anomalous magnetic moment, a¢??—a?M = (25.5+8.0) x 10~1° [L6], for which 
we require the SUSY effects to account at 20 level. (We note that 30 effects are considered to be inconclusive 


in high energy physics. In collider experiments, there are a large number of channels and observables and there 
is a good chance that some of the measurements can show such deviation from expectation. The muon g — 2 
experiment is quite different because there is just one quantity to measure in the experiment. In our opinion, 
the significance of the deviation should be taken rather seriously.) 


(6) Dark matter constraints from the WMAP relic density 0.0945 < Qh? < 0.1287 [17] and CDMS II limits on the 
scattering cross section [s]. To show the effects of the CDMS II limits, we will display the results with/without 
such limits. 


In addition to the above experimental limits, we also consider the constraint from the stability of the Higgs potential, 
which requires that the physical vacuum of the Higgs potential with non-vanishing vevs of Higgs scalars should be 
lower than any local minima. Further, the soft breaking parameters are required to be below 1 TeV to avoid the 
fine-tuning, and À (at weak scale) is less than about 0.7 to ensure perturbativity of the theory up to the grand 
unification scale (A_is increasing with the energy scale Eg). Note that most of these constraints have been encoded 
in NMSSMTools (1.4). We extend _this package and use it in our calculations. For the cross section of -nucleon 
scattering, we use the formulas in ig, for the MSSM and extend them to the NMSSM/nMSSM (see Appendix A). 
Considering all the constraints listed above, we scan over the parameters in the following ranges 


100 GeV < (MRY, MEP", ma, p) < 1 TeV, 
50 GeV < Mı < 1 TeV, 1< tang < 40, 
(Al lel) < 0.7, |Ax| < 1 TeV, (14) 


To reduce the number of the relevant soft parameters, we work in the so-called m°” scenario with following choice 
of the soft masses for the third generation squarks: Mo, = Mu, = Mp, = 800 GeV, and X; = A; — p cot 8 = —1600 
GeV. The advantage of such a choice is that other SUSY parameters are easy to survive the constraints (so that 
the bounds we obtain are conservative). Moreover, we assume the grand unification relation for the gaugino masses: 
Mı : Mə : M3 ~ 1 : 1.83 : 5.26 and also assume universal masses Mz and Mg for the three generations of sleptons and 
the first two generations of squarks respectively. 
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FIG. 1: The scatter plots for the spin-independent elastic cross section of y-nucleon scattering. The ‘+’ points (red) are 
excluded by CDMS limits (solid line), the ‘x’ (blue) would be further excluded by SuperCDMS 25kg in case of unobservation 
(dash-dotted line), and the ‘o’ (green) are beyond the SuperCDMS sensitivity. 


The surviving points are displayed in Fig. [I] for the spin-independent elastic cross section of .-nucleon scattering. 
We see that for each model the CDMS II limits can exclude a large part of the parameter space allowed by current 
collider constraints and the future SuperCDMS (25 kg) can cover the most part of the allowed parameter space. 
For the MSSM and NMSSM the dark matter mass range m,.o is from 50 GeV to 400 GeV, while for the nMSSM 


the dark matter mass is constrained below 40 GeV by current experiments and further constrained below 20GeV by 
SuperCDMS in case of unobservation. For the MSSM/NMSSM the LSP lower bound at 50 GeV is from the chargino 
lower bound of 103.5 GeV plus the assumed GUT relation Mı œ 0.5Mo; while the upper bound at 400 GeV is from 
the bino nature of the LSP (M; cannot be too large, must be much smaller than other relevant parameters) plus the 
constraints from the LEP II search for Higgs bosons, the muon g-2 and B-physics. Note that if we do not assume the 


GUT relation Mı ~ 0.5M2, then Mı can be as small as 40 GeV and the LSP lower bound in MSSM/NMSSM will 
not be sharply at 50 GeV. 
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FIG. 2: Same as Fig. [| but projected on the plane of [Nui /? and |Nis 2 versus dark matter mass. 


In Fig. 2] we show the bino component of x} in MSSM/NMSSM and the singlino component of x}? in nMSSM. We 
see that for both the MSSM and NMSSM x? is bino-dominant, while for the nMSSM x? is singlino-dominant, and the 
region allowed by CDMS limits (and SuperCDMS limits in case of unobservation) favor a more bino-like x? for the 
MSSM/NMSSM and a more singlino-like x? for the nMSSM. For the MSSM/NMSSM, the reason is obvious because 
the dominant contribution to the cross section comes from Fig. [9]in the Appendix and a more bino-like x? tends to 
suppress not only f in Eq. (83) [20], but also fË by diminishing Thoo. As for the nMSSM, xÌ is singlino-like due to 
the small singlino mass in the neutralino mass matrix. The peculiarity of the nMSSM predictions will be discussed 
at the end of this section. 


500 


400 


300 | 


u (GeV) j u (GeV) i u (GeV) 


FIG. 3: Same as Fig. [1| but projected on the plane of Mı versus u. The dashed lines are for Mı = p. 


In Fig. B]we project the surviving points on the plane of Mı versus u. We see that for both the MSSM and NMSSM 
most of the survived points are below the Mı = p line, implying that x? is bino-dominant. The region allowed by 
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the CDMS limits tends to have a larger u, indicating a more bino-like y?, which can be inferred from the neutralino 
mass matrices in Eq. (I). For the nMSSM the upper bound of 500 GeV for p is from the fact that a larger u leads 
to a lighter LSP (as will be shown in Eq{I5), which is then constrained by the required annihilation rate of the LSP. 

In Fig. [4] we display the surviving points on the plane of the chargino mass m + versus myo. For both the MSSM 
and NMSSM, the CDMS limits tend to favor a heavier chargino and ultimately the SuperCDMS limits tend to favor 
a wino- dominant chargino with mass about 2m,o. This can be understood because the CDMS/SuperCDMS limits 


require a large u, which makes yf to be dominated by wino with a mass M2 ~ 2M1. 
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FIG. 5: Same as Fig. [1| but projected on the plane of u versus the charged Higgs mass. 


In Figs] and [6] we display the surviving points on the plane of u and tan 8 versus the charged Higgs mass. In 
both the MSSM and NMSSM, large u and small tan 8 are favored for a light charged Higgs boson. The reason is 
as follows. In the MSSM, there are two CP-even Higgs bosons contributing to the cross section. One is the SM-like 
Higgs boson h? with mass around 120GeV and the other is the heavy boson H° with mass nearly degenerate with 
the charged Higgs boson. Then from the expression of f{7 in Eq.(A3), one can learn that the H° contribution to the 
scattering cross section get enhanced for light charged Higgs boson. In this case, to alleviate such enhancement, large 
u (to lower THoo) and/or small tan 8 (to lower Thq:q;) are needed. In the NMSSM, although there are three CP-even 
Higgs boson contributing to the scattering, we can get the same conclusion as the MSSM because one of the bosons 
is singlet-dominant and its contribution is suppressed by Thqiqi, and the contributions from the other two bosons are 
quite similar to the case of the MSSM, . 
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FIG. 6: Same as Fig. [1| but projected on the plane of tan 8 versus the charged Higgs mass. 
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FIG. 7: Same as Fig. [1] but projected on the plane of |A| versus the charged Higgs mass in NMSSM and nMSSM. 


In Fig. [7] we show the value of |A| versus the charged Higgs mass in NMSSM and nMSSM. This figure indicates 
that A larger than 0.4 is disfavored for the NMSSM. The underlying reason is that Thoo in Eq.(A3) depends on A 
explicitly and large À can enhance Thoo ug. By contrast, although CDMS has excluded some points with large A in 
the nMSSM, there are still many surviving points with A as large as 0.7. 

In Fig. [8] we show the decay branching ratio of h? + x9 versus the mass of the SM-like Higgs boson h°. Such a 
decay is strongly correlated to the y-nucleon scattering because the coupling h°y?9 is involved in both processes. We 
see that in the MSSM and NMSSM this decay mode can open only in a very narrow parameter space since x? cannot 
be so light, and in the allowed region this decay has a very small branching ratio (below 10%). By contrast, in the 
nMSSM this decay can open in a large part of the parameter space since the LSP can be very light, and its branching 
ratio can be quite large (over 80% or 90%). Such a large invisible decay ratio may indicate a severe challenge for 
finding the Higgs boson h? at the LHC if the nMSSM is the true story. Fig. B]also indicates that the mass of h? can 
reach 160 GeV. We checked that these cases correspond to A varying from 0.6 to 0.7 so that the mass is enhanced at 
tree level. 

Now we discuss the reason for the peculiarity of the nMSSM predictions shown from Fig. 2] to Fig. [8] About the 
narrow parameter space of the nMSSM constrained by collider experiments, a detailed analysis has been given in [4], 
here we only explain the behavior of the nMSSM under the CDMS/SuperCDMS limits. Our explanation is based 
on following three facts. The first comes from the neutralino mass matrix in Eq.(L1) which implies that myo can be 
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FIG. 8: Same as Fig. [I] but projected for the decay branching ratio of h° > x$ x? versus the mass of the Higgs boson h (the 
SM-like Higgs boson). 


written as [4]: 


2A? (v2 + v3) tan 8 
2u? + A? (v2 + v2) tan? B+ 1 


~ 
Mo = 


(15) 
This formula shows that to get a heavy x}, we need a large A, a small tan 8 as well as a moderate u. The second fact 
is that, due to the singlino dominance of x? in the nMSSM, the interaction of x? with squarks is suppressed and the 
Higgs mediated contribution in Figs[9]and[10]then becomes dominant in the scattering. In this case, A determines the 
size of the scattering for a given myo by affecting the coupling Thoo and a large À can enhance the cross section. 
The last fact is based on Fig. []which shows that the constraints of CDMS results become stringent for heavy X} and 
as a result, only myo around 40 GeV is excluded by CDMS. With these facts, one can easily understand the features 
of Figs. 28] For example, Fig. [6]and Fig. [/] indicate that the disfavored points by CDMS are characterized by small 
tan 8 and large A. The reason is that only under these two conditions, both mg and the cross section can be large 
simultaneously. 

The similarity of the allowed parameter space for the MSSM and NMSSM can be understood as follows. In both 
models x? is composed dominantly by bino, as shown in Fig. P} Then the properties of x9 (like the relic density and 
the_y-nucleon scattering) are similar in both models. Our such conclusion agrees with except that the conclusion 
of is based on a different scan scheme. Compared with foi, we considered more constraints and so our conclusions 
are more robust. 


IV. SUMMARY 


Considering the current direct and indirect collider constraints, we gave a comparative study for the neutralino dark 
matter scattering on nucleon in the MSSM, the NMSSM and the nMSSM. We showed the predictions for the elastic 
cross section by scanning over the parameter space allowed by the collider constraints and demonstrated the property 
of the allowed parameter space with/without the new limits from CDMS II. We found that for each model the new 
CDMS limits can exclude a large part of parameter space allowed by current collider constraints. The property of the 
allowed parameter space is found to be similar for MSSM and NMSSM, but quite different for n MSSM. Further, the 
future SuperCDMS can cover most part of the allowed parameter space for each model. 
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Appendix A: Spin-Independent Cross section of ,-nucleon scattering 
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FIG. 9: Feynman diagrams contributing to the scalar elastic-scattering amplitude of a neutralino from quarks in the MSSM, 
where H and h denote the CP-even Higgs bosons and ĝ represents a scalar quark. 


FIG. 10: Feynman diagrams contributing to the gluonic interaction with neutralinos, which contributes to the scalar elastic- 
scattering amplitude for neutralinos and nuclei. 


In supersymmetric models, the spin-independent elastic y-nucleon scattering is described by the following effective 
ue coin (20): 


2 
+basxxGt, Ge” — as(Bip + Bis)XO,0,x GOH” 
+asB2sX (ibp Ww a iy Yu) X gO, (A1) 


where the twist-two quark and gluon operators are defined by 


POSS AN 1 eats 
L = faXxdat Iq 2:7" ə” xO"), — sma AX | 


if. 7 1 | 
ol, =3 fonda + EWI — 59nd an ; 


1 
(2) uv a apv ald Vale 
Gee = GG a +99 G Gp 


10 


GH” is the gluon field-strength tensor, and f,, gg, b and B are coefficients. 
In the MSSM, the coefficients f, and qq are determined by calculating the diagrams in Fig. [9] in the extreme 
nonrelativistic limit, and they are given by 


5 1 Xa ijoWa ijo gThooThqiqi 
fa = AAi ewe Moon, 
7 g 4 > må, — (my + Mg)? 2 2m? 
2 2 
Ja. = SS (Xp igo) + (Weiso) (A3) 
i 7? 
a |m, — (my + ma)? 


where the subscripts q = u,d and i = 1,2,3 refers to the flavor index in quark sector, and X49, W4:;0; 
Thqiqı are the coupling coefficients of qi PrXoĝj, q PRXoĝj, XoXoh and qqh vertices respectively. 


The coefficients of the last four operators can be obtained in a similar way from Fig. [I0] and their expressions are 


Thoo and 
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m m. 
b = Tå Bp + Bs a ag ot Bis), 
1 g Thoo Gnas Gi 
De OG. ee) 
967 E n? m, g Tu 
1 
Bp = gap De aAa gijo Ii (Mq, Må; My), 
qi 45 
il 1 2 
Bs = 307 x5 ( 1 ago) + ( 740) In(™mq,,Mg;,Mx), 
qi sj 
1 
Bip = T57 2 Ma Xg ij 0WG igo Fa (Masa; My); 
qi 45 
1 1 2 2 
Bis = Ton L x5 K ago) + (Wy ago) | Ta(TMq,,™g,, mx), 
qi sj 
1 1 ! 2 ! 2 A 
Eng = ae 5 | (Xo) + (W; go) Is (Mai, Mg; Mx), aa 
qi 45 


where Ips are functions given in Eqs. (B.1a-e) of with the Eq. (B.1d) corrected as follows: the factor in the first 


term should read (m3 — m? — m2), with a corrected exponent for my; the term immediately following should read 


—1/ mm4, again with a corrected exponent for my; finally, a sign in the last term should be corrected so that it reads 


[e m3 + m2] L. 

About above formulae, two points should be noted. One is the Lagrangian in Eq. (AI) is specified at a high- 
energy scale, for example, Ho ~ mp, and in order to get the scattering rate measured in dark matter direct detection 
experiments, one must consider important QCD and SUSY-QCD corrections to the coefficients(2]}. In our calculation, 
we have considered such effect. The other is some extensions of the MSSM, such as NMSSM and nMSSM considered 
in this paper, usually predict extra CP-even Higgs bosons and neutralinos, and consequently, the couplings appeared 
in above formulae may be changed. In this case, the formulae listed above still keep valid in the sense that one must 
use the corresponding new couplings with the same convention as that in {19} and also include the contributions from 
new intermediate states. For example, the NMSSM predicts three CP-even Higgs bosons, and one should add the 
three boson contributions in getting fg (23). 

Given the effective Lagrangian in Eq.(Al), one can write down the spin-independent scattering cross section of a 
neutralino from a nucleon N (proton or neutron) in a standard way om) 


4m? 
oS! = A (A5) 


Mx MN 
mMmy +My 


where M, = is the reduced LSP mass, and fy is the effective couplings of the neutralino to nucleon, which 
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is given by: 


3 
~ 9x >D 9q(H0)4" (Ho) — Tè fFe 


+576" (u8)a45(18) [Bas + TZ (Bin + Bis)). (A6) 


In Eq.(A6), A denotes the fraction of the nucleon mass my that is due to the light quark q, and ffe = #(1 — 


JA, — fA — fA) is the heavy quark contribution to my, which is induced via gluon exchange. The function q^ (2) 
and G™ (u2) appear in the second moment of the quark (gluon) distribution functions and they represent the quark 
and gluon densities in the nucleon at the scale uo. The quantities gq(j10)q® (ug) and GY (u2)as(u2) in the third term 
and the last term is a renormalization-group invariant (in other words, independent of jo) and their evaluation was 
described in detail in [20]. In our calculation, we use ory = 64 MeV and o9 = 35 MeV to get the values of foa and 
use CTEQ6L to get the values of q (m?) and GN (mĉ). 

Before we end this section, we remind two subtleties in Eq. (AG) (20). One is to get the coefficient b by the formula 
in Eq.(A4), one should not include the contribution of u,d,s quarks to Bp since they are non-perturbative effects. 
The other is only top quark contribution needs to be considered in getting B2s in the last term. The reason is the 
contributions from u, d,s quarks to Bəs are non-perturbative effects, and the contributions from c, b quarks have been 
moved to the third term of Eq.(A6). 
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